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ABSTRACT 

We report the discovery of IIAT-P-27b, an exoplanet transiting the moderately bright G8 dwarf star 
GSC 0333-00351 {V = 12.214). The orbital period is 3.039586 ± 0.000012 d, the reference epoch of 
transit is 2455186.01879 ± 0.00054 (BJD), and the transit duration is 0.0705 ± 0.0019 d. The host star 
with its effective temperature 5300 ± 90 K is somewhat cooler than the Sun, and is more metal-rich 
with a metallicity of +0.29 ± 0.10. Its mass is 0.94 ± 0.04 Mq and radius is 0.90t^:g| Rq. For the 

planetary companion we determine a mass of 0.660 ± 0.033 Mj and radius of 1.038lQQ5g -^J- ^'-'^ 
30 known transiting exoplanets between 0.3 Mj and 0.8 Mj, a negative correlation between host star 
metallicity and planetary radius, and an additional dependence of planetary radius on equilibrium 
temperature are confirmed at a high level of statistical significance. 

Subject headings: planetary systems — stars: individual (HAT-P-27, GSC 0333-00351) — techniques: 
spectroscopic, photometric 



1. INTRODUCTION 

Studying exoplanets is vital for understanding our own 
Solar System, particularly its formation. The sample 
of more than 500 confirmed exoplanctf[3 so far enables 
us, for example, to test accretion and migratio n theo- 
ries ([Ida fc Lii]|[2008l ). study tidal interac tions tMardlinai 
120071) . examine atmospheric structures (jFortnevI I2O1O0 . 
and investigate correlations between the existence of 
planetary comp anions and the host star's metallicity 
(jida fc Linll2004l ) , and between the mass of close-in plan- 
ets a nd the spectral type of their host star (Ida & Lin 
l2005l) . 

Among these planets, transiting ones are of special sig- 
nificance, because the transit parameters yield planetary 
mass and radius estimates. They also provide a means 
to determine some of the stellar parameters more accu- 
rately than is possible with spectroscopy alone, such as 
the stellar surface gravity. The more than 100 transit- 
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ing exoplanets confirmed to date provide a sample large 
enough to draw meaningful conclusions about the plane- 
tary parameters that could not be determined by radial 
velocity (RV) data alone; for example, the correlation be- 
tween stellar chromospheric activity and planetary sur- 
face gravity (|Hartmanll2010( ) , or the correlation of plane- 
tary parameters with host star metallicity and planetary 
equilibrium temperature, as described in Section 21 

The Hu ngarian- made Automated Telescope Network 
(HATNet; iBakos et all 1201 If ) is a system of fully au- 
tomated wide-field small telescopes designed to detect 
the small photometric dips when exoplanets transit their 
host stars. Since 2006, HATNet has announced and pub- 
lished 26 planetary systems with 28 planets in total, 26 
of which transit their host stars. 

Here we report the detection of our twenty-seventh 
transiting exoplanet, named HAT-P-27b, around the rel- 
atively bright G8 dwarf known as GSC 0333-00351. This 
planet is a textbook example of a transiting exoplanet 
with its radius 1.04 i?j and orbital period 3.0396 d being 
close to the median values for currently known transiting 
exoplanets, and with its mass of 0.66 Mj being not much 
less than the median mass of transiting exoplanets. 

In Section [2] we present the photometric detection of 
the transit, along with photometric and spectroscopic 
follow-up observations of the host star HAT-P-27. In 
Section [31 we describe the analysis of the data, first ruling 
out false positive scenarios, then determining parameters 
of the host star, and finally performing a global fit for all 
observational data. We conclude the paper by discussing 
HAT-P-27b in the context of other known transiting ex- 
oplanets and investigating correlations of planetary pa- 
rameters with host star metallicity and equilibrium tem- 
perature in Section 131 

2. OBSERVATIONS 
2.1. Photometric detection 
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Transits of HAT-P-27b were detected in two HAT- 
Net fields containing its host star GSC 0333-00351, 
also known as 2MASS 14510418+0556505; a = 
14'^51'^04.32" 5 = -h05°56'50.5", J2000, V=12.214 
proege et alj 120061: hereafter HAT-P-27. These fields 
were observed in Sloan r-band on a nightly basis, weather 
conditions permitting, from 2009 January to August, 
with the HAT-6 and HAT-10 instruments on Mount Hop- 
kins, and with the HAT-9 instrument on Mauna Kea. In 
total, we took 10 600 science frames with 5 minute ex- 
posure time and 5.5 minute cadence. For approximately 
1200 of the images, photometric measurements of indi- 
vidual stars had significant error, therefore these frames 
were rejected. Each image contains approximately 20 000 
stars down to r w 14.5. For the brightest stars, the per- 
point photometric precision is approximately 4.5 mmag. 



dates and the results of these observations are summa- 
rized in Tabl e [H 
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Fig. 1. — Top panel: Unbinncd photometric data of HAT- 
P-27 consisting of 9 400 Sloan r-band 5.5 minute cadence mea- 
surements obtained with HATNet telescopes, folded with period 
P = 3.039586 d. A simple transit curv e fit to the data points 
is displayed with a solid line. See Section 13.31 for details. Bottom 
panel: A close-up view of the transit. Small gray dots are the same 
folded data as above; large black dots show the light curve binned 
in phase using a bin size of 0.002. 

Calibration, astrometry, aperture photometry, Exter- 
nal Parameter Decorrelation (EPD), the Trend Filter- 
ing Algorithm (TEA) and the Box-fitting Least Squares 
method we re appl ied to the data as described in 
iBakos et all (|2010aD . We detected a transit signature 
in the light curve of HAT-P-27, with a signal depth of 
10.6 mmag and period of P = 3.0396 days. This pre- 
sumed transit has relative first-to-last-contact duration 
q — 0.0232 ± 0.0006, corresponding to total duration 
0.0705 ± 0.0019 days (1.693 ± 0.046 hours). The folded 
light curve is presented in Figure [TJ 

2.2. Reconnaissance Spectroscopy 

Reconnaissance spectra (Lath am et al.l l2009a') of HAT- 
P-27 were taken using three facilities : the Tillingh ast Re- 
flector Echelle Spectrograph (TRES: 11X1^120081) on the 
1.5 m Tillinghast Reflector at FLWO, the echelle spec- 
trograph on the 2.5 m du Pont telescope at Las Cam- 
panas Observatory (LCO) in Chile, and the echelle spec- 
trograph on the Australian National University (ANU) 
2.3 m telescope at Siding Spring Observatory (SSO) in 
Australia. We gathered two spectra of HAT-P-27 with 
TRES in 2009 July and 2010 February, two spectra with 
the du Pont telescope in 2009 July, and fourteen spectra 
with the ANU 2.3 m telescope in 2009 July. The exact 



Following lOuinn et all (|2010D and iBuchhave et al.l 
(201C|), we calculated the mean radial velocity, effec- 
tive temperature, surface gravity, and projected rota- 
tional velocity of the host star, based on spectra taken 
by TRES. The inferred radial velocity RMS residual of 
0.05 kms""'^ is consistent with no detectable RV variation 
within the precision of the measurements. We estab- 
lished the following atmospheric parameters for HAT-P- 
27: effective temperature Toff* = 5250 ± 100 K, surface 
gravity logg* = 4.75 ± 0.25 (cgs), and projected rota- 
tional velocity wsini = 1.5 ± 1.0 kms~^, indicating a 
G8 dwarf. The mean heliocentric RV of HAT-P-27 af- 
ter subtracting the gravitational redshift of the Sun is 
= -15.765 ± 0.51 kms'^. 
Because this is the first time we used the du Pont 2.5 
m and ANU 2.3 m telescopes for reconnaissance spec- 
troscopy of HATNet targets, we briefly describe the in- 
struments and our data reduction procedure. 

The spectrograph on the du Pont telescope was used 
with a 4" long and 1.5" wide slit. The obtained spectra 
have wavelength coverage « 3700-7000 A at a resolution 
of A/AA « 26 000. During the first observation the see- 
ing ranged between 2-3" and we used an exposure time 
of 1200 s, which provided ~ 3000 electrons per resolu- 
tion clement at the wavelength of 5187 A. The seeing 
during the second observation was « 1.8" and we used 
an exposure time of 150 s to obtain a lower S/N spectrum 
sufficient to detect a velocity variation of several kms^^. 
We obtained a ThAr lamp spectrum before and after 
each observation to use in determining the wavelength 
solution. We used the CCDPROC package of IRAIof to 
perform overscan correction and fiat-fielding of the im- 
ages, and the ECHELLE package to extract the spectra 
and to determine and apply the dispersion corrections. 

The extracted du Pont spectra were then cross- 
correlated against a library of synthetic stellar spectra to 
estimate the effective temperature, surface gravity, pro- 
jected rotational velocity, and radial velocity of the star. 
We followed a pro cedure similar to that described by 
[Torres et all (pOOl . using the same synthetic templates, 
but broadened to the resolution of the du Pont echelle. 
These templates, which were ge nerated for the CfA Digi- 
tal Speedometer (|Lathamlll992[ ) , only cover a wavelength 
rang e of 5150-5360 A, so we restricted our analysis to a 
single order of the spectrum covering a similar range. 

Spectra were also taken using the echelle spectrograph 
on the ANU 2.3 m telescope. The echelle was used 
in a standard configuration with a 1.8" wide slit and 
300 mm~^ cross-disperser setting of 5°50', which deliv- 
ered 27 full orders between 3900-6720 A with a nominal 
spectral resolution of A/AA « 23 000. The CCD is a 
2Kx2K format with 13.5 /im x 13.5 /im pixels. The 
gain was two electrons per ADU resulting in a read noise 
of approximately 2.3 ADU for each pixel. The spectra 
were binned by two along the spatial direction. A to- 
tal of fourteen 1200 s exposures were taken. The seeing 
ranged from 2" to 3". The signal-to- noise on a single 

IRAF is distributed by the National Optical Astronomical 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 
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TABLE 1 

Summary of reconnaissance spectroscopy observations of HAT-P-27 



Instrument 


Date 


Number of 






V sin i 


7rv^ 






Spectra 


(K) 


(cgs) 


(kms-i) 


(kms-i) 


TRES 


2009 Jul 05 


1 


5250 


4.5 


2 


-15.75 


du Pont 


2009 Jul 10 


1 


5250 


4.0 





-16.03 


du Pont 


2009 Jul 11 


1 


5000 


3.5 





-18.03 


ANU 2.3 m 


2009 Jul 16 


5 








-20.58 


ANU 2.3 m 


2009 Jul 17 


6 








-21.32 


ANU 2.3 m 


2009 Jul 18 


3 








-20.61 


TRES 


2010 Feb 13 


1 


5250 


5.0 


1 


-15.78 



^ The mean heliocentric RV of the target. Systematic differences between the veloci- 
ties from different instruments arc consistent with the velocity zero-point uncertainties. 



pixel was between 5 and 10 for each individual exposure. 
ThAr lamp calibration exposures were taken every hour 
for wavelength calibration. A high signal-to-noise expo- 
sure was also taken of the bright radial velocity standard 
star HD 223311. 

Spectra were reduced using tasks in the IRAF pack- 
ages CCDPROC and ECHELLE. The spectra were 
cross-correlated against the radial velocity standard star 
HD 223311 using the IRAF task FXCOR in the RV pack- 
age. We used at least 20, typically 25 of the 27 orders for 
the cross-correlation, rejecting the bluest orders for many 
exposures where the signal-to-noise was too low. Each 
spectral order was treated separately and the mean of 
the velocities from the individual orders was calculated. 
Their standard deviations were less than 0.65 kms^^ for 
each exposure. 

Each night, the exposures were taken within a two 
hour interval, much shorter than the orbital period of 
the presumed companion. For detecting large radial ve- 
locity variations to rule out the possibility of an eclips- 
ing binary, we consider the mean radial velocities per 
night. The standard deviations between exposures were 
less than 0.75 kms~^ for each night. Stellar parameters 
could have been estimated only with large uncertainty 
based on data with such low signal-to-noise. Therefore 
these parameters are not calculated from the ANU 2.3 m 
observations. 

The results of the observations taken with these three 
telescopes are listed in Table [1] Note that for each tele- 
scope, the radial velocity measurement uncertainty is 
much higher than the radial velocity variations of the 
Sun due to Solar System bodies. Therefore we only cal- 
culated heliocentric radial velocities of the target. For 
the more accurate measurements described in the next 
section, we will use barycentric radial velocities instead. 
This accuracy, however, is enough to rule out the case of 
an eclipsing binary star with great certainty. The largest 
RV variation within an instrument was only 2 kms~^, 
much less than the orbital speed in a typical binary sys- 
tem. Note that the zero-point shift between instruments 
is as large as '--^ 5 kms~^, due to the different methods 
used for data reduction. Also, all eighteen spectra were 
single-lined and spectral lines were symmetric, providing 
no evidence for additional stars in the system up to the 
precision of the measurements. 

2.3. High resolution, high S/N spectroscopy 

We acquired high-resolution, high-S/N spectra of 
HAT -P-27 using the HIRES instrument (jVogt et al.l 
|1994() on the Keck I telescope located on Mauna Kea, 



Hawaii, between 2009 December and 2010 June. The 
spectrometer was configured with the 0.86" wide slit, 
yielding a resolving power of A/AA « 55 000 over the 
wavelength range of « 3800-8000 A. 

Nine exposures wer e taken using an I2 gas cell (see 
iMarcv fc Bu"tle^ [1991 . and a single template exposure 
was obtained witho ut the absorption cell. We followed 
iButler et all (|1996[) to establish RVs in the Solar Sys- 
tem barycentric frame. We a lso calculated the S index 
for each spectrum (following llsaacson" fc Fisched [20101. 
The resulting values and their uncertainties are listed in 
Tabic [21 They are plotted period- folded in Figure [21 to- 
gether with the fit established in Section [3l 

The effective temperature esta blished later in 
Sectio n 13.21 compared to Figure 4 of IValenti fc Fischeil 
(|2005f) implies B - V = 0.800 for th e star. This can in 
turn be used in the formula given in iNoves et al.l ()1984[ ) 
together with the median ^hk of 0.231 to conclude 



logi?^ 



HK 



-4.785. This activity value is consistent with 



the RV jitter value 6.3 ms establishe d in Sect ion 13.31 
according to the obser vations given bv I Wright! ^005.) . 
Based on Figure 9 in llsaacson fc Fisched ( 2010f ). this 



value qualifies HAT-P-27 as chromospherically active 
relative to California Planet Search targets of the same 
spectral class. The activity index does not correlate 
significantly with orbital phase. 

2.4. Photometric follow-up observations 

A high-precision photometric follow-up of a complete 
transit was carried out, permitting refined estimates of 
the light curve parameters and thus orbital and planetary 
properties: The transit of HAT-P-27 was observed on the 
night of 2010 March 2 with the KeplerCam CCD camera 
on the FLWO 1.2 m telescope. We acquired 168 science 
frames in Sloan i-band with 60 second exposure time, 73 
second cadence. 

Following the procedure described by iBakos et al.l 
pOlOail . these images were first calibrated, then astrom- 
etry and aperture photometry was performed to arrive 
at light curves, which were finally cleaned of trends us- 
ing EPD and TEA, carried out simultaneously with the 
global modeling described in Section 15751 The result is 
shown in Figure [3l along with the best-fit transit light 
curve; the individual measurements are reported in Ta- 
ble [3 

3. ANALYSIS 
3.1. Excluding blend scenarios 
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TABLE 2 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-27 



BJD (UTC) RV^' anv^ BS'^ ubs S'^ Phase 

(2400000+) (ms-i) (ms-i) (ms^i) (ms-^) 



55192.13748 






-4.96 


3.92 


0.211 


0.013 


55193.14273 


-87.86 


1.73 


-2.89 


5.30 


0.212 


0.344 


55194.14692 


73.31 


1.51 


-4.20 


2.47 


0.204 


0.674 


55252.00750 


86.53 


1.75 


4.37 


2.03 


0.227 


0.710 


55257.15592 


-51.27 


1.84 


1.11 


3.44 


0.244 


0.404 


55261.02101 


84.94 


2.02 


-0.91 


3.63 


0.247 


0.675 


55290.06295 


-84.71 


1.89 


3.14 


3.57 


0.245 


0.230 


55312.92780 


98.63 


2.08 


-6.69 


7.87 


0.252 


0.752 


55374.90154 


-74.42 


1.85 


-10.88 


19.39 


0.230 


0.141 


55375.82176 


-29.80 


2.10 


-3.76 


7.04 


0.233 


0.444 



Note. — For the iodine-free template exposures there is no RV measure- 
ment, but the BS and S index can still be determined. 



^ The zero-point of these velocities is arbitrary. An overall offset 710! fitted 
to these velocities in Section l3.3l has not been subtracted. 

^ Internal errors excluding the component of astrophysical/instrumcntal jitter 
considered in Section l3.3l 

^ The bisector spans have been corrected for sky contamination following 
IHartman et all 120111) . 

^ Relative chromospheric activity index, calibrated to the scale of 
IVaughan, Preston fe WilsonI I I1978I '). 



TABLE 3 

High-precision differential photometry of HAT-P-27 



BJD (UTC) Mag'' o-Mag Mag(orig)'' Filter 
(2 400 000-I-) 



55258, 
55258, 
55258, 
55258, 
55258, 
55258, 
55258, 
55258, 
55258, 
55258, 



88089 
88157 
88242 
88327 
88412 
88499 
88584 
88689 
88773 
88859 



0.00295 
0.00160 
-0.00094 
-0.00140 
0.00104 
0.00093 
-0.00041 
0.00056 
0.00123 
0.00098 



0.00109 
0.00109 
0.00095 
0.00094 
0.00095 
0.00094 
0.00095 
0.00095 
0.00094 
0.00094 



10.86800 
10.86810 
10.86430 
10.86320 
10.86640 
10.86500 
10.86470 
10.86350 
10.86580 
10.86540 



Note. — This table is available in a machine-readable form 
in the online journal. A portion is shown here for guidance 
regarding its form and content. 

^ The out-of-transit level has been subtracted. These mag- 
nitudes have been obtained by the EPD and TFA procedures, 
carried out simultaneously with the transit fit. 
^ Raw magnitude values without application of the EPD and 
TFA procedures. 



To further exclude possible blend s, we perform bisec - 
tor analysis the same way as in §5 of lBakos et al.l (|2007| ). 
A significant scatter is found, strongly correlated to the 
presence of moonligh t, which we account f or using the 
method described by IHartman et al.] ()2011[ ). The bisec- 
tor spans, corrected for the effect of the moonlight, are 
shown in the third panel of Figure [51 They do not exhibit 
significant correlation with the RV values, and the RMS 
scatter of the bisector spans (4.6 ms~^) is significantly 
smaller than the RV amplitude. These findings rule out 
a blend scenario with high certainty, implying that the 
measured photometric and spectroscopic features are due 
to a planet orbiting HAT-P-27. 

3.2. Properties of the parent star 

We first determine spectroscopic parameters of HAT- 
P-27, which will allow us to calculate stellar mass and 
radius. The Spectroscopy Made Easy analysis pack- 



age (SME; iValenti fc Piskunovl I1996D is used to estab- 
lish the effective temperature, metallicity and stellar sur- 
face gravity based on the Keck/HIRES template spec- 
trum of H AT-P-27, using atomic line data from the 
database of lValenti fc Fischeii (|2005D . After an initial es- 
timate for these parameters, we perform a Monte Carlo 
calculation, relying also on the normalized semi-major 
axis a/ Ri, inferre d from transi t light c urves, for the rea- 
sons described bv lBakos et al.l (|2010bf ). The final values 
adopted after two iterations are T^s* = 5300 ± 90 K, 
[Fe/H] = -1-0.29 ± 0.10, and vsini = 0.4 ± 0.4 kms-\ 
also listed in Table 21 

Based on the final spectroscopic parameters the model 
isochrones yield a stellar mass M^, = 0.945 ± 0.035 AIq 
and radius = 0.898t|^:^^9 ^r HAT-P-27, along 
with other properties listed at the bottom of Table IH 
These values classify the star as a G8 dwarf, and sug- 
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Fig. 2. — Top panel: Keck/HIRES RV measurements for 
HAT-P-27 shown as a function of orbital phase, along with our 
best-fit model (see Table (5)1. Zero phase corresponds to the time of 
mid-transit. The center-of-mass velocity has been subtracted. Sec- 
ond panel: Velocity O—C residuals from the best fit. The error bars 
include a component from astrophysical/instrumental jitter (6.3 
ms~^) added in quadrature to the formal errors (see Section l3.3p . 
Third panel: Bisector spans (BS), with the mean value subtracted, 
and corrected for sky contamination. The measu rement from the 
template spectrum is included (see Section l3.f|l . Bottom panel: 
Chromospheric activity index S. Again, the measurement from 
the template spectrum is included. Note: Panels have different 
vertical scales. The data point replotted in the second period is 
represented by an open symbol. 
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Fig. 3. — Unbinned Sloan j-band transit light curve, acquired 
with KeplerCam on the FLWO f .2 m telescope on 2010 March 2. 
The light curve has been EPD and TFA-processed. Our bes t fit 
from the global modeling is shown by the solid line. See Section 13.31 
for details. Residuals from the fit are displayed at the bottom. The 
error bars represent the photon and background shot noise, plus 
the readout noise. 



lYi et all ()2001[ ) for a metallicity of +0.29 are pfotted in 
Figure m along with the best estimate of a/i?^ and T^e* 
of HAT-P-27, and their la and 2a confidence ellipsoids. 
For comparison, the initial SME result, corresponding to 
a somewhat younger state, is also indicated. 




5500 
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Fig. 4. — Model isochrones from IYi et al.l 1 120011 ) for the measured 
metallicity of HAT-P-27, [Fe/H]= +0.29, and ages running from 
1 Gyr to 14 Gyr in 1 Gyr increments, left to right. The adopted 
values of T^g-i, and a/Ri, are shown as the solid dot, surrounded 
by Icr and 2cr confidence ellipsoids. The initial values of Tpfj^ and 
a/Rt, from the first SME and light curve analyses are represented 
with a triangle. 

The intrinsic absolute magnitude predictions of this 
model (given in the ESO photometric system) can 
be compared to observations to calculate the dis- 
tance of HAT-P-27. For this we use the near-infrared 
brightness measureme nts from the 2MASS Catalogue 
(jSkrutskie et alj [20Q6f ): J2MASS = 10.626 ± 0.026, 
-ff2MASs = 10.249 ±0.023 and JtT^MAss = 10.10 9 ± 0.021. 
These values are conver ted to ESO (|Carpen tcr 2001), 
then compared to the absolute magnitude estimates to 
calculate the distance. We account for interstellar dust 
extinction in the line of sight using E{B~V) = 0.036 ± 
0.010 from the dust map bv lSchlegel elaP ([T991 FI with 
a conservative uncertainty estimate. This has to be mul- 
tiplied by a factor depending on the distance of the star 
and its Galactic latitude (see iBonifacio et "aI][200Q ). As- 
suming diffuse interstellar medium and no dense clouds 
along the line of sight, we use the value R y = 3.1. along 
with t he coefficients given in Table 3 in iCardelli et al.l 
(|1989f l. These let us calculate extinction parameters for 
each band based on the reddening. Finally, comparing 
extinctions, absolute magnitude predictions and 2MASS 
apparent magnitudes for J, H and K bands, we arrive 
at separate distance estimates. These are in good agree- 
ment, yielding an average distance of 204 ± 14 pc. The 
uncertainty does not account for possible systematics of 
the stellar evolution model. Note that this value is only 
1 pc less than the more simple estimate ignoring extinc- 
tion. The model predicts an unreddened color index of 
{J —K)i-aodc\ — 0.50±0.04. Reddening would change it to 
an estimated observed value of (J — if)red = 0.514 ± 0.04 
using the above parameters. This matches the actual 
observed color index (J — i^)ESO = 0.550 ± 0.037 within 
la. 

3.3. Global modeling of the data 



gest an age of 4.4^2^ Gyr. The model isochrones of 
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TABLE 4 
Stellar parameters for HAT-P-27 



TABLE 5 
Orbital and planetary parameters 



Parameter Value Source 

Spectroscopic properties 

Tcff* (K) 5300 ± 90 SME'' 

[Fe/H] +0.29 ±0.10 SME 

i;sini (kms-i) . . . 0.4 ± 0.4 SME 

iJmac (kms-i) . . . . 3.29 SME 

Vmic (kms-i) .... 0.85 SME 

7RV (kms-i) -15.765±0.51 TRES 

Photometric properties 

V (mag) 12.214 TASS 

V-Zc (mag) 0.527 ±0.12 TASS 

J (mag) 10.626 ± 0.026 2MASS 

H (mag) 10.249 ± 0.023 2MASS 

(mag) 10.109 ±0.021 2MASS 

Derived properties 

(Mq) 0.945 ±0.035 YY+a/R*+SME 

R„ (Rq) 0.898i;iJ-g^^ YY+a/iJ^+SME 

logs* (cgs) 4.51 ±0.04 YY+a/i?*+SME 

(Lq) 0.57l°'o? YY+a/R*+SME 

Mv (mag) 5.55 ±0.17 YY+a/R*+SME 

Mx (mag, ESO) . 3.62 ± 0.12 YY+a/R*+SME 

Age (Gyr) 4.4l^-^ YY+a/iJ^+SME 

Distance (pc) 204 ± 14 YY+a/i?*+SME 

^ SME — Spectroscopy Made Easy package for the analysis of high- 
resolution spectra (Valenti & Piskunov 1996). These parameters rely 
primarily on SME, but have a small dependence also on the iterative 
analysis incorporating the isochrone search and global modeling of 
the data, as described in the text. 

^ YY+a/_R*+SME = Based on the YY isochrones JYi et al.llSoOll '). 
a/R^ as a luminosity indicator, and the SME results. 

We fi t the combined model described bv iBakos et al.l 
(|2010aD to the HATNet photometry, follow-up photom- 
etry, and radial velocity measurements simultaneously. 
The eight main parameters describing the model are the 
time of the first observed transit center, the time of the 
follow-up transit center, the normalized planetary ra- 
dius p = Rp/Ri,, the square of the impact parameter 6^, 
the reciprocal of the half duration of the transit C/^*: 
the RV semi-amplitude K , and the Lagrangian elements 
k = e cos uj and h = e sin uj (where uj is the longitude of 
periastron). 

Instrumental parameters of the model include the 
HATNet out-of-transit magnitude and the relative zero- 
point of the Keck RVs. The joint fit provides the full a 
posteriori probability distributions of all variables (in- 
cluding loggi,), which are used to update the limb- 
darkening coefficients for another iteration of the joint 
fit. This leads to estimated distributions for the stellar, 
light curve, and RV parameters, which are combined to 
calculate values for planetary parameters. These final 
values are summarized in Table [S] 

The orbital eccentricity is consiste nt with zero (using 
the method of lLucv fc Sweenevll 19711 we find that there is 
a 25% conditional probability of detecting an eccentricity 
of at least 0.078 given a circular orbit and an uncertainty 
of 0.047). Nevertheless, we stress that throughout the 
global modeling, the orbit was allowed to be eccentric, 
and all system parameters and their respective errors in- 
herently contain the uncertainty arising from the floating 
k and h values. 

4. DISCUSSION 



Parameter Value 

Light curve parameters 

P (days) 3.039586 ± 0.000012 

Tc^' (BJD, UTC) 2455186.01879 ± 0.00054 

Ti4 (days) ^ 0.0705 ± 0.0019 

Ti2 = T34 (days) 0.0214 ± 0.0053 

<^/R* 9.65t° « 

C/iJ* 38.20 ±0.80 

flp/fl* 0.1186 ±0.0031 

0.690in^^ 

b = acosi/Ri, 0.83llo;o2o 

i (degree) 84.7io * 

Limb-darkening coefficients 

ai (linear term) 0.3627 

bi (quadratic term) 0.2816 

RV parameters 

K (ms-i) 96.1 ±4.5 

fcRv'= 0.036 ±0.031 

/iRv'= 0.066 ±0.048 

e 0.078 ±0.047 

u) (degree) 63 ± 64 

RV jitter (ms"!) ' 6.3 

Secondary eclipse parameters 

Ts (BJD) 2455187.608 ± 0.060 

14 0.0739 ± 0.0061 

Ts'i2 0.0368 ± 0.0066 

Planetary parameters 

Mp (Mj) 0.660 ±0.033 

C(Mp,i?p) g 0.310 

Pp (gcm-3) 0.73 ± 0.13 

\oggp (cgs) 3.18 ± 0.05 

a (AU) 0.0403 ± 0.0005 

Tcq (K) 1207 ± 41 

0.054 ± 0.004 

Fper (108 ergs-1 cm-2) i . . . S.OllJ'j* 

Fap (10** ergs-i cm-2) i ^ ^ ^ 4 gg ± o.48 

{F) (10« ergs-i cm-2) i ^ ^ ^ ^ 4.79^°;^^ 

^ Reference epoch of mid transit that minimizes the correlation 
with the orbital period. 

^ Total transit duration, time between first to last contact. 

Ingress/egress time, time between first and second, or third and 
fourth contact. 

Adopted from the tabulations bv I Claret l|2004r according to the 
spectroscopic (SME) parameters listed in Table [4] 

Lagrangian orbital parameters derived from the global modeling, 
and primarily determined by the RV data. 

^ The contribution of the intrinsic stellar jitter and possible instru- 
mental errors that needs to be added in quadrature to the calculated 
RV uncertainties so that x^/dof — 1 in the joint fit. 
^ Correlation coefficient between the planetary mass and radius. 
^ The Safronov number is given by iHansen fc BarmanI 120071) as 
e = i(Vesc/V<„b)^ = (a/flp)(Mp/M,) . 

^ Stellar irradiation flux per unit surface area at periastron, apas- 
tron and time-averaged over the orbit, respectively. 

4.1. Properties of HAT- P- 27b 

Figure [5] presents the currently known transiting exo- 
planets and Solar System gas planets on a mass — radius 
diagram, with HAT-P-27b highlighted. Also s hown are 
the planetary isochrones of lFortnev et al.l (|2007[ ) interpo- 
lated to the insolation of HAT-P-27 at the orbit of HAT- 
P-27b. Taking into consideration the age established in 
Section 13.21 the planetary parameters are consistent with 
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a hot Jupiter with a 10 core in a 4 Gyr old system. 

HAT-P-27b can be seen to lie inside the large accu- 
mulation of planets with similar masses and radii. To 
further compare it to other Hot Jupiters, in Figure |6] we 
present histograms of mass, radius and period for the 112 
transiting exoplanets confirmed to date. 

When comparing these parameters, we note that there 
is only one transiting exoplanet known that is more mas- 
sive and has a smaller radius and a smaller period than 
HAT-P-27b: this is HAT-P-20b with 7.246 Mj, 0.867 i?j 
on a 2.875 d orbit (Bakos ct al. 2010b). This means that 
HAT-P-27b is less inflated than other planets of similar 
mass and orbital period, possibly due to a larger than 
average core. 

Regarding eccentricity, there are 31 transiting exoplan- 
ets known under 8 Mj with an orbital period within 0.5 
d of that of HAT-P-27b, out of which 8 - more than a 
quarter of them - arc thought to be eccentric. This hints 
that there is a possibility for the orbit of HAT-P-27b to 
be eccentric as well, justifying our choice not to fix ec- 
centricity to zero in Section [S] Future observations of 
radial velocity or occultation timing would be required 
to determine whether the orbit is indeed eccentric. 

The im pact parame t er of HAT-P-27b is unusually 
large. As iRibas et al.l ()2008D pointed out, such a near 
grazing transit has the advantage of its depth and dura- 
tion being more sensitive to the presense of further plan- 
etary companions on inclined orbits. This makes HAT- 
P-27 a promising target for transit timing variation and 
tr ansit duration v ariati on studies. 

iKnutson et al.l ()2010[ ) found a strong negative corre- 
lation between chromospheric activity of the host star 
and temperature inversion in the planetary atmosphere. 
However, since early type stars dominate magnitude lim- 
ited surveys, cool, that is, active planetary hosts are rare. 
The bottom right panel of Figure [6] shows that HAT- 
P-27 is relatively active compared to known planetary 
hosts for which log i?HK been reported, making it an 
exciting target for Spitzer Space Telescope to test this 
correlation. 

4.2. Correlation of planetary parameters with host star 
metallicity 

The relation between host star metallicity (denoted as 
[Fe/H]^, for clarity, not to be confused with the assumed 
metal content of the planet) and planetary composition 
was studied by Guillot ct al. (2006). A positive correla- 
tion, with Pearson correlation coefficient r = 0.78, was 
found between the inferred mass of the planetary core 
and stellar metallicity for the seven transiting exoplanets 
known at that time with positive inferred core mass. The 
idea is that planets have formed from the same cloud as 
their host stars, and therefore their metal content should 
correlate. However, it is not clear how stellar metallicity 
is connected to planetary metallicity, especially because 
a larger rocky core is likely to accrete more gas during 
t he planet's formation. 

iBurrows et al.l (|2007l ) also investigated this relation, 
based on 12 transiting exoplanets known at the time. 
They used an atmospheric opacity dependent core mass 
model to explain radius anomalies. They also found a 
strong correlation between host star metallicity and in- 
ferred core mass, but the correlation coefficient was not 
reported. 



lEnoch et al.l (|2010D found that there is a strong neg- 
ative correlation with r = —0.53 between [Fe/H]^, and 
Rp for the 18 known transiting exoplanets below 0.6 
Mj, whereas this correlation is negligible for more mas- 
sive planets. This can be explain ed by notici ng that 
the theoretical pla n et mo dels of iFortnev et a .1 (12007') . 
iBodenheimer et all (|2003D . and iBaraffe et al.l (|2008[r an 
suggest that the radius of a planet is more sensitive to 
its composition for low mass planets than it is for more 
massive ones. 

In this subsection, we examine further the correlation 
between host star metallicity and planetary mass or plan- 
etary radius. We use a substantially expanded sample of 
30 known transiting exoplanets with masses between 0.3 
Mj and 0.8 Mj, see Table ID The upper limit is selected 
to exclude more massive planets whose radius is expected 
to depend less on the composition, see above. We explain 
the role of the lower limit and the effect of the two lowest 
mass planets in Table [H] later. 

The null hypothesis is that the host star metallicity and 
the selected planetar parameter are independent. The al- 
ternative hypothesis is that they are related by some un- 
derlying phenomenon. A false positive, also known as an 
error of the first kind, is rejection of the null hypothesis in 
spite of it being true. We implement three independent 
statistical methods to estimate the false positive proba- 
bility: t-test, bootstrap technique and F-test. We denote 
the probability estimates by pi, p2 and pa, respectively. 
This is the statistical significance of the correlation: the 
lower this probability is, the more confidently the null 
hypothesis (i.e., no correlation) can be rejected. 

For the i-test, we assume that the investigated parame- 
ters have normal distribution. For each set of data pairs, 
we calculate the t value from the correlation coefficient r 
and sample size n using the formula 



t 



1 — 



The conditional distribution of this variable given the 
null hypothesis is Stude nt's t distrib ution with n — 2 de- 
grees of freedom (.Press et al.l 11993 , p. 640). Then the 
estimate pi for false positive probability is determined 
by looking up the two-tailed probability of this distribu- 
tion yielding larger absolute value that the one measured. 
For comparison, we also performed the t-test for the sam- 
ples and par a meter s studied by Guillot ct al. ( 2006) and 
lEnoch et al.l (|2010D . The resulting values are listed in 
Table [71 

For the sample set of Table [SI we also implement 
the b ootstrap resampling technique ([Efron fc Tibshiranil 
I2003D . This has the advantage that no assumption 
about the a priori distribution of the data is neces- 
sary. To perform bootstrap resampling, consider the data 
(xi,yi), (x2,?/2), • ■ • , (a;„,2/„), where is the host star 
metallicity, and yi is the corresponding planetary mass 
or radius. Again, we would like to calculate an estimate 
P2 of the probability that a sample of similar distribution, 
but independent parameters for each pair, has a correla- 
tion coefficient that exceeds that of our measurements in 
absolute value. For this, we build 10 000000 sample sets 
of n pairs by drawing x and y values independently with 
replacements from the set of measured x and y values, 
respectively. The percentile rank of the absolute value of 
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TABLE 6 

Parameters of 30 transiting exoplanets between 0.3 Afj and 0.8 Afj in increasing order of mass 



Mp{Mj) 



Rp {Rj ) 



7-off(K) 

1262 ± 31 
17921^3^ 

1565t|J 

1417+62 
1628lg3 

1438 ± 39 
1568 ± 33 
1020 ± 17 
16621112 
1194^5^ 
1306 ± 30 

787 ± 15 
1652 ± 28 
1025"* 



[Fe/H], 



reference 



WASP-21b 
HD 149026b 

Kepler-7b 

WASP-13b 
Kepler-8b 

CoRoT-5b 
WASP-31b 

WASP-ll/HAT-P-lOb 
WASP- 17b 
WASP-6b 
HAT- P- lb 

HAT-P-17b 
WASP-15b 
OGLE-TR-lllb 



HAT-P-4b 

WASP-22b 
XO-2b 

HAT-P-25b 
WASP-25b 
WASP-34b 
HAT-P-3b 

HAT-P-28b 
HAT-P-27b 
HAT-P-24b 
HD 209458b 

Kepler-6b 

OGLE-TR-lOb 
CoRoT-4b 
TrES-lb 
HAT-P-9b 



0.3 ±0.01 
0.368lO:»li 



1.07 ±0.05 

n O1O-I-0.027 



0.416 



+0.036 
0.035 



1.439 



+0.058 
0.056 



r.+0.056 
-0.05 



0.463 

0.46 ±0.14 



1.21 
1.31 



+0.14 
-0.12 
+0.076 
0.08 



n 467+°-''*^ 
U.4D/_Q Q24 

0.478 ± 0.03 

0.487 ±0.018 
40+0.059 
'^•*^-0.056 

0.524 ±0.031 

0.53 ±0.019 
0.542 ±0.05 
0.55 ±0.1 



1 Qoo+0.046 

1.000_Q 

1.54 ±0.06 

i.UUO_Q 

1.74t«i« 
1.224t0;«^^ 
1.225 ±0.059 



1293 ± 0.03 



1.428 
1.019 



+0.077 
0.077 

+0.026 
0.026 



0.556 ±0.068 1.367^0.044 



0.56 ±0.02 

-+0.055 
-0.055 



0.566"' 



1.12 ±0.04 

j+0.029 
-0.028 



0.983"* 



0.567 ±0.022 

0.58 ±0.04 

0.59 ±0.01 

596"*"°"2* 
u.oyo_o 026 

0.636 ± 0.037 
0.660 ± 0.033 
0.685 ± 0.033 



1.19^ 



0.685 



+0.015 
-0.014 



K0.081 
-0.056 

1 r,r, + 0.06 

-^•^^-0.05 
1 99+0.11 
-^•^^-0.08 

899''"° ''*^ 
U.033_Q 049 

i.l»y_0 075 
i.U3»_o 058 

1.242 ±0.067 

1.0d3_Q 019 



0.62 



+0.025 
0.028 



1.164 



+0.025 
■0.017 



0.62 ±0.14 
0.72 ±0.08 

7t;2+0-0*7 
U. /OZ_Q 045 

0.780 ± 0.090 



-+0.14 
J-0.12 
F0.06 
-0.05 
7+0.022 
-0.021 

1.40 ±0.06 



1.25: 
1.191 
1.067"* 



h26 
-25 



+30 



1686 



1430 ± 30 
1319i^* 



-0.4 ±0.1 
+0.36 ±0.05 

+0.11 ±0.03 

+0.0 + 0.2 
-0.055 + 0.03 

-0.25 + 0.06 
-0.19 + 0.09 
+0.13 + 0.08 
-0.25 + 0.09 
-0.20 + 0.09 
+0.21 + 0.03 

+0.0 + 0.08 
-0.17 + 0.11 
+0.19 + 0.07 



+0.24 + 0.08 



-0.05 + 0.08 
0.44 + 0.04 



Bouchv et al. (2010) 
Ammler-vonEiffetaL 
Carter et al. (2009) 
Latham et al. (2009b), 
Kipping fc Bakos (20101 
Skillen et al. (20091 
Jenkins et al. (20101 . 
Kipping fc Bakos (20101 
Rauer et al. (2009) 
Anderson ct al. (20101 
Bakos et al. (2009) 
Anderson et al. (2010 1 
Gillon et al. (20091 
Torres et al. (2008), 
Ammlcr-von Eiff e t al. 
Howard ct al. (2010T~ 




Torres et al. (2008), 
Winn ct al. (2010) 
Maxted ct al. (2010) 
Torres et al. (2008), 
Ammler-von Eiff ct al. 



1202 + 36 


+0, 


.31 


+ 


0, 


.08 


Ouinn et al. (20101 


1212 + 35 


-0, 


.07 


+ 


0, 


.1 


Enoch et al. (20101 


1250 + 30 


-0, 


.02 


+ 


0, 


.1 


Smallev et al. (20101 


1127"*"*'-* 
iJ-^' -39 


+0, 


.27 


+ 


0, 


.04 


Torres et al. (20071. 












Torres et al. (20081 


1371 + 50 


+0, 


.12 


+ 


0, 


.08 


Buchhave et al. (20111 


1207 + 41 


+0, 


.29 


+ 


0, 


.1 


this paper 


1637 + 42 


-0, 


.16 


+ 


0, 


.08 


Kipping et al. (20101 


1449 + 12 


+0, 


.01 


+ 


0, 


.03 


Laughlin et al. (2005), 


1459^24 












Torres ct al. (2008)" 


+0, 


.34 


+ 


0, 


.04 


Dunham ct al. (201^), 












Kipping & Bakos (2010 


1481±55 


+0, 


.15 


+ 


0, 


.15 


Torres et al. (20081 


1074 + 19 


+0, 


.05 


+ 


0, 


.07 


Moutou et al. (20081 


11401" 


+0, 


.02 


+ 


0, 


.05 


Torres et al. (20081 


1530 + 40 


+0, 


.12 


+ 


0, 


.2 


Shporer et al. (20091 



the correlation coefficient for the measured data among 
the random samples gives our estimates P2 , listed in Ta- 
ble H 

Finally, we test these cor relations with an additional 
method, the F-test (see e.g. lLuptonlll993l p. 100). This 
requires that the null hypothesis (no correlation) be 
nested in the tested hypothesis (linear correlation) , which 
indeed is the case. Let RSSi denote the residual sum of 
squares for the best fit of the null hypothesis, that is, the 
variance of yi about its mean, and RSS2 denote the resid- 
ual sum of squares of the linear fit. The no correlation 
model has = 1 free parameters: the mean, whereas the 
linear fit has I2 = 2. The conditional distribution of 

RSS1-RSS2 

~ RSS2 

n—l2 

given the null hypothesis is the F-distribution with {I2 — 
h, n—h) degrees of freedom. This enables us to calculate 



P3, the third estimate for the false positive probability. 

The estimates Pi, P2, P3 given by the three methods are 
listed in Table [T] For each correlation, they coincide up 
to the uncertainty of the methods. The values reflect the 
significa nt [Fe/H]^ — M n nrB an d [Fe /H]^ — Rr, correlatio ns 
found bv lGuillot et~an pOOl and lEnoch etall (pOlOl ). 

As for the 30 planets listed in TablelHl it is important to 
note that the correlations depend strongly on the choice 
of the lower mass limit. The two least massive plan- 
ets in the table are WASP-21b with a mass of 0.3 Mj, 
very low host star metallicity of —0.4, and average radius 
of 1.07 Rj; and the dense HD 149026b with a mass of 
0.368 Mj, high host star metallicity of +0.36, and low 
radius of 0.813 Increasing the lower mass limit for 
our sample first excludes WASP-21b, which would much 
support the positive [Fe/H]^, — Mp correlation with its 
low mass and low host star metallicity. Further increas- 
ing the limit then excludes HD 149026b, which would 
much weaken it with its low mass and high host star 
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metallicity. To have an unbiased result, outliers cannot 
be excluded without a justified reason, therefore we need 
to compare the false positive probabilities of the three 
nested samples. They scatter between 15% and 58%, 
neither supporting, nor rejecting a [Fe/H]^ — Mp correla- 
tion. 

Similarily, both WASP-21b and HD 149026b have a 
strong effect on the negative [Fe/H]^ — Rp correlation, 
because of the extreme value of their host star metallic- 
ities. In this case, we see that the maximum of the false 
positive probabilities is 0.44%, therefore this correlation 
is statistically significant for all our choices of lower mass 
limits. This is at least a fivefold im p rovem ent over the 
sample investigated bv lEnoch et al.l (|2010[ ). due to the 
larger sample size. 

4.3. Dependence on planetary equilibrium temperature 
Other factors like insolation are li kely to influ- 



ence pl anetary radiu s as w e ll, see e.g. iFqrtnev et al.l 
J2007), 'Kovac set all ([Moh . lEnoch et al.l ()2010t ). and 
[fiicdi ct al. (20"TT|). To further investigate this relation, 
we compare two models: for null hypothesis, we accept 
the linear planetary radius — host star metallicity relation 
of the previous section: 



[Fe/H],. 



(1) 



The second model - alternative hypothesis - is sim- 
ilar to that of lEnoch et al.l (|2010( ) , accounting for the 
equilibrium temperature Teq in addition to the host star 
metallicity: 



[Fe/H]. 



(2) 



The equilibrium temperature of the planet is calculated 
from the time-averaged stellar flux on its orbit, assuming 
gray body spectrum for the planets, and neglecting tidal 
and other heating mechanisms. For simplicity, we now 
include all 30 planets of Table [H] in our models. With the 
best fit parameters, the two models are 

^p = 1.235 i?j - 0.478 i?j • [Fe/H]^, 



0.690 i?j - 0.431 i?j • [Fe/H]^ + 0.398 Rj 



T, 



cq 



1000 K 



To quantify the statistical significance, we apply the 
F-test again. The resulting false positive probability 
is 0.18%. This means that once we accept the depen- 
dence of planetary radius on host star metallicity, then 
the probability of such a correlation with planetary equi- 
librium temperature if they were not physically related is 
only 0.18%. This strongly supports the three-parameter 
linear fit model. For reference, the residual sums of 
squares for the one, two and three-parameter fits of the 
[Fe/H]*— i?p data are 1.13 Rj, 0.84 Rj and 0.58 Rj, re- 
spectively. 



-0.505 between Rp and metal- 



Figure ?? presents Rp — /3^^ ■ T^q versus metallicity for 
the 30 planets. Equation ([2]) predicts this quantity to 
be R^^ + a" • [Fe/H]*, which is also plotted. HAT-P-27b 
apparently follows the model's prediction. For reference, 
the correlation coefficient between the displayed trans- 
formed variables is now r^^ = —0.536, which has a larger 
absolute value than 
licity, as expected. 

This analysis supports the statement that planetary 
radius depends on equilibrium temperat ure in a d dition 
to h ost star metallicity , as found by En och et al.l (|20100 
and iFaedi et al.l (|2011[ ). However, this correlation does 
not imply that insolation itself would inflate planets: the 
underlying phenomenon could be related to anything cor- 
related to equilibrium tempera ture, or eq uivalcntly, or- 
bital radius. For instance, Bat vein fc Ste venson (201(|) 
suggest that it is Ohmic dissipation in the interior of the 
planet that i nflates hot Jupiters. T his theory is further 
supported bv lLaughlin et al.l ()201lD . 

Altogether, HAT-P-27b is an important addition to 
the growing sample of low-mass Jupiters. It orbits a 
metal rich star, and supports the suggested correlations 
between host star metallicity, planetary equilibrium tem- 
perature, and planetary radius. Also, HAT-P-27 is chro- 
mospherically active, providing an excellent case for re- 
fining the confidence level of the hypothesized correla- 
tion between stellar activity and planetary temperature 
inversion. 
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Fig. 7. — Planetary radius corrected for linear equilibrium temperature dependence in Jupiter radii versus host star metallicity for the 
30 known transiting cxoplancts with masses between 0.3 M,j and 0.8 Mj, including HAT-P-27b (solid black circle), other HATNet planets 
(solid gray circles), and other known transiting exoplanets (solid gray squares, WASP-21b and HD 149026b labeled); with best linear fit 
overlaid (dashed line). 



